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Abstract: We study the phenomenology of mixed-sneutrino dark matter in the Minimal 
i?-symmetric Supersymmetric Standard Model (MRSSM). Mixed sneutrinos fit naturally 
within the MRSSM, as the smallness (or absence) of neutrino Yukawa couplings singles out 
sneutrino ^-terms as the only ones not automatically forbidden by /^-symmetry. We per- 
form a study of randomly generated sneutrino mass matrices and find that (i) the measured 
value of r^DM is well within the range of typical values obtained for the relic abundance 
of the lightest sneutrino, (ii) with small lepton- number- violating mass terms rn^^hh for 
the right-handed sneutrinos, random matrices satisfying the il-DM constraint have a decent 
probability of satisfying direct detection constraints, and much of the remaining parameter 
space will be probed by upcoming experiments, (Hi) the rn^^fin terms radiatively generate 
appropriately small Majorana neutrino masses, with neutrino oscillation data favoring a 
mostly sterile lightest sneutrino with a dominantly ///r-fiavored active component, and {iv) 
a sneutrino LSP with a significant /U component can lead to striking signals of e-/i flavor 
violation in dilepton invariant-mass distributions at the LHC. 



1. Introduction 

The weak scale promises to hold a wealth of new physics. Both the hierarchy problem and 
the success of WIMP dark matter candidates suggest that new particles may be discovered 
at the Tevatron or LHC. Indeed, one of the appealing features of supersymmetric theories 
with R-parity is that the lightest superpartner (LSP) is a natural dark matter candidate. 

At the same time, the minimal supersymmetric standard model (MSSM) is fraught 
with issues that invite us to explore new and interesting structures. In particular, the 
flavor problem of the MSSM is significant [1, 2], requiring a suppression of flavor violating 
terms at the 10~^ level, even with CP conservation. 

The Minimal i?-symmetric Supersymmetric Standard Model (MRSSM) [3] was pro- 
posed as a new solution to the supersymmetric flavor problem without flavor-blind me- 
diation. By imposing a continuous or suitably large discrete i2-symmetry on weak scale 
supersymmetry, order one flavor violating soft masses for the sfermions are allowed. The 
i?-symmetry forbids Majorana gaugino masses, trilinear A-terms and the ;U-term. Since 
massless gauginos and Higgsinos are in conflict with experiment, the MSSM is augmented 
by considering /^-symmetry preserving Dirac gauginos whose masses are generated by "su- 
persoft" operators [4]. Although previous attempts to implement Dirac gauginos [5, 6, 7] 
were within the context of flavor-blind SUSY breaking, in the MRSSM they are considered 
part of a general softly broken supersymmetric theory. 

The phenomenology of the MRSSM is striking. Because the flavor violation is order 
one, dramatic flavor- violating signatures [8, 9, 10, 11] are possible. The TeV-scale Dirac 
gauginos and 0(1) flavor violation in the sfermion sector open up the possibility of a natural 
mixed sneutrino dark matter candidate [12, 13] as well as interesting flavor violating signals 
at the LHC [14]. 

The MSSM sneutrino was long ago considered an intriguing dark matter candidate 
[15, 16] but is today ruled out by the combined relic abundance, direct-detection, and 
invisible Z-width constraints. By mixing the sneutrino with a gauge-singlet field, the 
coupling of the lightest mass eigenstate to the Z can be suppressed and the invisible Z- 
width constraint is met for mixing angles satisifying sin0 < 0.4 [12, 17, 18], for arbitrarily 
small masses. The appropriate relic abundance can also be achieved since the sneutrino 
annihilation rate in the early universe is suppressed [12]. Related scenarios of sneutrino 
dark matter are explored in Refs. [19, 20, 21, 22, 23]. 

However, mixed-sneutrino dark matter is still tightly constrained by direct detection 
experiments[12], even with the suppression of the Z-coupling. Splitting the scalar and 
pseudo-scalar components can relax direct-detection constraints [24, 25, 12], but once ex- 
perimental upper bounds on the radiatively generated neutrino masses are taken into ac- 
count, the viable parameter space is rather finely tuned [23, 13]. This makes it challenging 
to realize a model of sneutrino inelastic dark matter (iDM) [26] in standard supersymmetric 
scenarios. 
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These tensions are resolved in the MRSSM, where mixed-sneutrino dark matter more- 
over finds a natural home. Because the flavor violation is naturally large, there is a natural 
understanding of the origin of the large neutrino mixing angles through an anarchic [27] 
setup. Because the gauginos are Dirac, with only small Majorana masses coming from the 
conformal anomaly, the radiative masses are naturally a loop factor lower than in simple 
extensions of the MSSM, and are thus appropriately small, even with large enough lepton- 
violation to evade direct detection constraints. Finally, since neutrinos have small (or zero) 
Yukawas, they are singled out as the only superfields in the MRSSM to have large A-terms, 
which are essential for radiatively generating the neutrino mass matrix. 

A study of mixed sneutrino dark matter in the MRSSM is thus strongly motivated. In 
this paper, we will explore the degree to which the anarchic scenario, implemented here, 
can produce neutrino oscillation parameters consistent with experimental results, while 
simultaneously satisfying direct detection constraints and yielding the appropriate relic 
abundance. The connection to neutrino masses will give us insight into the likely flavor 
content of the sneutrino LSP, which will in turn motivate analyses sensitive to lepton flavor 
violation at the LHC. Because the radiative corrections to neutrino masses are suppressed 
due to the Dirac gauginos, mixed-sneutrino iDM becomes viable again. If the scenario we 
explore is correct, we may learn that the one form of non-baryonic dark matter that we 
already know (neutrinos) is intimately linked to the other (cold dark matter). 

2. Mixed sneutrinos in the MRSSM 

In addition to the elimination of ^-terms, the /i-term and Majorana gaugino masses, the 
MRSSM has additional fields that distinguish it from the MSSM. Adjoint chiral super- 
fields are added to marry the standard gauginos. The Higgs sector is extended by adding 
multiplets Ru and Rd with the appropriate charges to allow i?-symmetric ;U-terms to be 
written down with Hu and Hd respectively. The scalar components of the Higgs (and 
not the i2-fields) acquire vevs that break electroweak symmetry, thereby preserving the 
i?- symmetry. 

Because the gauginos are Dirac, they can naturally be much heavier than the scalars 
[4]. The supersymmetric flavor problems are addressed by the TeV-scale Dirac gauginos 
and absence of j4-terms at small tan /3 and by the extended Higgs sector at large tan/3. 
Flavor violating observables in AF = 2 processes are suppressed due to two effects. Ra- 
diative corrections to squark masses from Dirac gauginos are a finite one-loop effect, which 
allows them to be naturally lighter than in the MSSM, and the heavy gauginos provide 
a suppression by a factor m?/m? ~ 10~^. The i?-symmetry allows leading flavor violat- 
ing operators [1, 2, 28] to be dimension 6 instead of 5 in the case of Majorana gauginos, 
thereby leading to another suppression m|/m| ~ 10"^. AF = 1 processes such as /i — )• 67 
or 6 — )• 57, involve a helicity flip in the diagram. For Dirac gauginos, the opposite helicity 
state has no direct couplings to the sfermions and these diagrams are absent [29]. Only 
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the much smaUer diagram with the hehcity flip on the external fermion line or internal 
line from Higgsino-gaugino mixing can be present. In the MSSM, the radiatively gener- 
ated couplings of down-type quarks to Hu at large tan /3 give the largest contributions 
to FCNCs [30]. However, these contributions are absent in the MRSSM due to different 
diagrams being eliminated by the absence of ^-terms, the /i-term and Majorana gauginos. 

Mixed sneutrino models arise naturally within the MRSSM. Since the holomorphic 
structure of the Yukawa and j4-terms is the same, i?-symmetry allows eii/ier Yukawa terms 
or A-terms to be written down, but not both. The electrically charged fermions of the 
Standard Model have large Yukawa couplings, and thus no A-terms are possible. In the 
neutrino sector, the Yukawas are zero, and thus large A-terms are natural. The trilinear 
interaction Aijhiljhu couples the left-handed sneutrinos to one or more additional singlet 
superfields Ni. The large flavor violation in the slepton sector leads to large flavor mixing 
among the left-handed sneutrinos. By including three singlet superfields Ni, we obtain a 
6x6 sneutrino mass matrix with maximal flavor mixing and large off-diagonal A-terms 
that mix the left and right handed sneutrinos. 

The inclusion of the A terms, restricting to one generation, leads to the mass matrix 

m| + cos 2/3 At; sin/3\ (2 1) 

Av sin /3 I 

where m^, is the soft scalar mass for the left-handed sleptons, ^rn^ cos 213 is the D-term 
contribution to the left-handed sneutrino mass-squared and m/j is the soft scalar mass for 
the gauge-singlet right-handed sneutrino. 

Although mi,m/j and A are independent parameters, it is natural to have m/j < 
in the low-energy theory [12]. Since n is a standard-model singlet, there are no gaugino 
loops to drive its mass upward as the energy scale runs down from Adpi to My[/. There are 
new, sizeable loop diagrams arising from the A-terms [Aijlifijhu + AijDihjhu) which push 
the soft masses down. Since two states, / and P, propagate in the loop contributing to mfi 
and only one, h, propagates in the loop contributing to m^, mn is driven down faster than 
mi- 

2.1 Radiative neutrino mass 

The AL=2 Majorana neutrino mass term, HuHuLlLl is i?-symmetric and allowed in 
the MRSSM [3]. In the presence of a lepton-number- violating mass term m^„nn for the 
right-handed sneutrinos, radiative neutrino mass can be generated via neutralino/sneutrino 
loops [31, 17, 32] 1. These dia grams require a Majorana mass insertion and are absent if 
the gauginos are Dirac. In the MRSSM, anomaly mediation naturally generates a small 
Majorana gaugino mass in addition to the large, 0{TeV) Dirac mass. Furthermore, since 

^To generate a m^„nn term we need _R-symmetry breaking to feed into the right-handed sneutrino 
masses. This gives us reason to expect this mass to be small - in the present framework we want it to be 
of order 100 MeV. 
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the Dirac gauginos are somewhat heavier than in the MSSM, the radiative neutrino mass 
is suppressed and allows small neutrino masses without tuning. 

The radiatively generated neutrino mass can be calculated in the mass insertion ap- 
proximation [33] . Choosing a basis for the fermions and bosons where the gauge couplings 
to the neutralinos are flavor diagonal, the flavor changing is exhibited by non-diagonality 
of the sfermion propagators. In the mass-eigenstate basis, the neutrino mass diagram has 
odd numbers of insertions from the Majorana neutralino mass and from Ajj, the lepton 
number violating mass term in the mass basis. The Majorana neutralino mass can be 
treated as a perturbation on the Dirac neutralino mass matrix. 

.„.. X ..-.„ 



> ^^/^^.AAAAA/WW^^vWW^ ' 

Ve Xl 

Figure 1: The radiatively generated Majorana neutrino mass in the mass eigenstate basis 




For a single wino Majorana mass insertion {"^y/ is the Dirac wino mass), the 
diagram is calculated as (Fig. 1): 



where 
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327r2 (m?- — m?)(m^- — m^)(m^ — m^) 

-,2 2 i„„r 2 /_2i , 2 



log["^vi//"^i ] + log["^i l^w^ + "^i log["^i h^j\ ) ' (2-3) 

9ei n ^L^ii „ __o9i, 



Ztt^ P'^ — — mf ^ p'^ — 

= ^^^^ ^ 1 + 2 2 log "^vk/"^^ ' 2.4 

^ W ^ ' \ w ^ / 

and i,j run over the 6 mixed sneutrinos. 

The neutrino mass matrix thus obtained can be diagonalized with a matrix U . Exper- 
imental results on neutrino masses and mixings can thus give insight into the structure of 
the sneutrino mass matrix [27]. 
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In our analysis below, we will consider experimental constraints on the following quan- 
tities from the neutrino sector: 

R = Am22/Ami3, (2.5) 

SC=4|C/e3|'(l-|C/e3|'), (2.6) 
Satm=4|C/^3|2(l-|C/^3|2), (2.7) 

Sq = A\Ue2?\Uel\\ (2.8) 

where Amfj is the smallest splitting and Am23 is the next largest splitting of the neutrino 
masses. 

The current experimental bounds on these quantities (at 3 a) are [34]: 

0.027 <R< 0.04, (2.9) 
sc < 0.05, (2.10) 
0.34 < Satm < 0.67, (2.11) 
0.26 < Sq < 0.4. (2.12) 

2.2 Direct detection 

Direct detection experiments constrain the spin-independent and spin-dependent WIMP- 
nucleon cross sections. In the presence of lepton-number violation, the mixed LSP sneutrino 
{vi) can evade direct detection constraints arising from from Z boson exchange [26, 35, 25]. 
However, the h-vi-vi coupling is modified due to the additional gauge degrees of freedom 
and also due to the large A-terms. Therefore, scattering from Higgs exchange still con- 
strains the theory [12]. In the decoupling limit {mA ^ "^h), the cross section for Higgs- 
mediated i>i-nucleon scattering is 

ao = ( ( ^^^y ( (2.48 X 10-^W), (2.13) 
Vl-26xl0-3y \mN + mi,J \ ruh J ^ ^ ' 

where ttt-tv is the nucleon mass, we take ghNN = 1.26 x 10~^ for the Higgs-nucleon coupling, 
with the up, down, strange and heavy quark contributions taken from [36, 37, 38], and 
9huiui is the coupling of the light Higgs boson to the LSP sneutrino Di. Recent lattice 
measurements of the light and strange quark condensates in the proton affect the cross- 
section evaluated here by ~ 20% at most [39]. 

The h-vi-vi coupling receives contributions from the Higgs L'-term and the large A- 
terms. For pure Dirac gaugino masses, the i^-term is zero and the Higgs quartic vanishes. 
Electroweak symmetry breaking can be recovered by adding small Majorana masses that 
lift the D-flat direction. The resultant coupling for the LSP sneutrino in the mass basis 
becomes 

^ ^ 2 2 \ 

2 , 2 1 0772 - 1 



(2.14) ^ 
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where m^i. and rriDi are the Majorana and Dirac gaugino masses, and U is the mixing 
matrix that relates the sneutrino mass (i) and gauge (a) eigenstates by = Uiai^a- In 
the formula above, k runs over the active flavors (e,/i,T) and p runs over the three sterile 
flavors of right-handed sneutrinos. 

Currently, the strongest direct detection bounds on WIMP-nucleon scattering are from 
CDMS [40] and XENONIO [41]. CDMS sets an upper limit on the spin-independent cross 
section of 4.6xl0~^^ cm^ at the 90% confidence level for a WIMP mass of 60 GeV/c^ and 
XENONIO has an upper bound of 4.5 x 10"^'' cm^ for a WIMP mass of 30 GeV/c^. 

2.3 Invisible 2^-viridth 

The invisible Z-width is extremely constrained. From current experimental bounds, T{Z — )• 
invisible) = 499.0 ± 1.5 MeV and r(Z vv) = 501.65 ± 0.11 MeV. Thus, new light degrees 
of freedom with a total mass less than Mz/2 that couple to the Z boson must have a decay 
rate within AT < 3 MeV [42]. 

The rate for Z decaying to two mixed sneutrinos is given by 

where U is the mixing matrix that relates the mass and gauge sneutrino eigenstates, k runs 
over e,/U,T and T^, = 167 MeV. The invisible Z-width constraint becomes 

^r,,,^, <Ar (2.16) 

where rrii + rrij < Mz- 

2.4 Relic abundance 

The dominant annihilation channels for mixed sneutrinos in the early universe are s-channel 
Z-exchange, t-channel neutralino exchange and s-channel Higgs exchange. The Higgs- 
sneutrino coupling is enhanced due to large j4-terms and the mixing suppresses the Z- 
exchange contribution so the s-channel Higgs exchange is often dominant. Sneutrino masses 
near the Higgs or Z poles should give a depleted relic abundance due to rapid annihilation. 

In the case of 6 mixed sneutrinos, the LSP sneutrino mainly contributes to the dark 
matter relic abundance. The heavier sneutrinos annihilate into the LSP sneutrino before 
freezeout. Co-annihilation is possible if the mass difference between two sneutrinos is small. 
The current bound on the dark matter energy density is 0.094 < QoMh'^ < 0.128 [43]. 

3. Analysis 

The sneutrino mass matrix in the MRSSM can have order one fiavor mixing in the left 
and right handed sectors as well as large off-diagonal A-terms. The large fiavor violation 
is shielded from and uncorrelated with the fiavor observables. The fiavor violating terms 
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and A-terms in the sneutrino mass matrix can be essentially random without contributing 
to the supersymmetric flavor problem. 

However, the sneutrino mass matrix is not unconstrained, due to limits from the ra- 
diatively generated neutrino masses and the invisible Z-width. As a plausible dark matter 
candidate, the sneutrino mass matrix is further constrained and must satisfy direct detec- 
tion and relic abundance bounds, as discussed above. 

We perform a systematic study of randomly generated 6x6 sneutrino mass matrices 
and analyze the likely properties of a mixed sneutrino dark matter candidate consistent 
with the experimental constraints discussed above. Our analysis sheds light on the mass 
hierarchies of the sneutrinos, whether the dark matter is likely to be predominantly active 
or predominantly sterile, and the expected flavor composition of the active part. 

We generate 75 million random 6x6 Hermitian (real, symmetric) sneutrino mass 
matrices with positive eigenvalues. Each sneutrino matrix has 6 random entries for m/, 
and mji, 9 for the ^-terms and 3 for A, the lepton number violating mass terms for the right 
handed sneutrinos. The matrix entries are randomly sampled from a uniform distribution. 
The entries for mi,mR are sampled from [0,250] GeV and the range for A is [—100, 100] 
GeV. We take the elements of A to be small perturbations on the sneutrino mass matrix, 
also with linear distributions. In our analysis the only effect of A is to generate neutrino 
masses, and only the ratios of the elements of A to one another are important. This is 
because the neutrino masses are also directly proportional to the Majorana gaugino masses, 
which we also treat as a small perturbations. Two cases are considered: Run I, with heavy 
gauginos and an additional loop-generated flavor diagonal mass of 200 GeV for the left 
handed sneutrinos (arising from gauge mediation), and Run II with light gauginos and 
without the flavor diagonal mass. 

To calculate the relic abundance of the mixed sneutrino LSP, we use the MicrOmegas 
2.0.7 code [44]. The randomly generated sneutrino matrix is related to a random slepton 
matrix by the left-handed slepton and sneutrino mass difference arising from the Z)-term 
contribution. Thus the left-handed sleptons are also mixed in this analysis. The relevant 
couplings of the lightest 4 mixed sneutrinos and the 3 left-handed sleptons are modifled 
by the appropriate mixing angles in the MSSM flies in CalcHEP [45] and fed into Mi- 
crOmegas to calculate the relic abundance. Since the lightest sneutrinos contribute to the 
relic abundance calculation, only the lightest 4 mixed sneutrinos are added to the code. 

Although the right-handed slepton mass matrix is also potentially mixed, we take it to 
be proportional to the identity matrix in this analysis. The right-handed slepton couplings 
do not affect the relic abundance calculation significantly, and restricting flavor violation 
to the left-handed sleptons leads to a clearer analysis of the collider signals discussed in 
Section 5. 

The gauginos are taken to have ©(TeV) masses. As we have seen in equation 2.14, 
the h-Di-Di contains free parameters due to the modified Higgs D-term. However, the 
contributions to this coupling from the ^-terms are numerically more important, and so 
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we can safely neglect the D-teim contribution to the coupling. 

We can now present the results of imposing the various experimental constraints on 
the sneutrino mass matrices. 

4. Results 

A table of the statistics obtained from applying neutrino mass, relic abundance, invisible 
Z- width and direct detection constraints is given below (Table 1). We generate 75 million 
sets of matrices, first for the case with a flavor-diagonal mass of 200 GeV (Run I - Heavy 
gauginos), and then for the case without a flavor-diagonal mass (Run II - Light gauginos). 

One interesting question is whether what we know about neutrino masses and mixings 
makes the scenario more or less viable. To answer this question, we first calculate the 
efficiency with which matrices passing the neutrino mass cuts subsequently satisfy the relic 
abundance, invisible Z-width and direct detection constraints. We then compare this to 
the efficiency with which random sets of matrices - ones not subjected to the neutrino mass 
cuts - satisfy the same remaining constraints. The "Random I" and "Random 11" rows 
of Table 1 show the results for 10,000 sets of random matrices unconstrained by neutrino 
measurements, for comparison with Run I and Run II, respectively. 

The results in Table 1 show that matrices which satisfy the neutrino mass constraints 
are about as likely to satisfy the remaining constraints as purely random matrices are. 
Thus, fitting to neutrino mass properties does not seem to benefit or impede the search. 



Run 










% events 










+r^_>.inv+direct 




I 


7754 


389 


343 


69 


0.89 


Random I 


10000 


428 


369 


106 


1.06 


II 


5757 


273 


144 


12 


0.21 


Random II 


10000 


309 


200 


27 


0.27 



Table 1: Results of experimental cuts on random sneutrino mass matrices, as explained in the 
text. 

4.1 Neutrino Masses and Mixings 

About one in 10^ of the random sets of of sneutrino mass matrices satisfy the neutrino mass 
constraints. Of course, this small efficiency is not by itself a discouraging result, because 
as the measurement errors associated with the various experimental constraints become 
smaller and smaller, the fraction of random matrices passing all cuts will inevitably tend 
toward zero. The question, then, is how should we evaluate the efficiency which which our 
matrices pass the neutrino mass cuts? 

Probably the simplest approach is to think of the cuts as tunings. That is, sc < 0.05 
is a 5% tuning, while the constraints on sq, Satm-, and R are 14%, 33%, and 1.3% tunings. 
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Figure 2: The sc distributions peak towards zero for the samples that satisfy the neutrino cuts, 
(i) Run I (ii) Run II 



Major ana 


no cuts 


Satm 


Sq 


Satm + S0 


no cuts 


75,000,000 


14,845,534 


8,069,129 


1,388,589 


sc 


9,007,887 


1,633,398 


760,902 


142,528 


R 


3,786,775 


405,762 


717,334 


67,116 


sc + R 


490,465 


84,996 


43,169 


7,754 



Table 2: The impact of neutrino mass and mixing cuts for Run I 



Majorana 


no cuts 


Satm 


Sq 


Satm ~l~ Sq 


no cuts 


75,000,000 


15,508,313 


7,861,926 


1,445,386 


sc 


9,488,149 


1,750,349 


816,193 


159,378 


R 


2,690,087 


284,677 


541,726 


51,241 


sc + R 


353,406 


62,412 


31,165 


5,757 



Table 3: The impact of neutrino mass and mixing cuts for Run II 



With this interpretation, we would expect only a fraction ~ 3 X 10~^ of the total to pass 
the cuts. The 10~^ which pass is (marginally) better than this, suggesting that there is 
nothing particularly surprising about the degree to which neutrino measurements constrain 
the sneutrino mass matrices. 

That said, the fact that a particular region of parameter space has been singled out 
allows us to make more specific statements about the nature of the LSP, as we will see 
below. 

Before we get to the LSP, let us consider the neutrino mass matrix in more detail. 
The radiative Majorana neutrino mass matrix generated via neutralino/sneutrino loops is 
also random. A table of the statistics obtained for the various neutrino mass and mixing 
cuts is given above (Tables 2 and 3). The neutrino mass matrices that pass the neutrino 
cuts generally exhibit a normal mass hierarchy, with less than 0.5% exhibiting an inverted 
hierarchy. 
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100 GeV 


sin 9 


0.2 


K 


1 TeV 





150 GeV 


6 


100 keV 


M 





Table 4: Benchmark point gives an experimentally viable neutrino mass, rrii, — 0.05 eV 

To check that the neutrino masses obtained are in the experimentally viable range, we 
consider a set of benchmark values for a 2 x 2 sneutrino matrix (Table 4) . sin 9 is the mixing 
angle and 6 is the mass splitting between the lightest CP-even and CP-odd sneutrino. A 
splitting of order 100 keV ensures that inelastic scattering at XENONIO is kinematically 
impossible. The Dirac Wino mass is heavy and the Majorana Wino mass is taken to be a 
one-loop factor down from the Dirac mass. In this case, the neutrino mass is given by [13] 



1 



327r2 



sin^ 95my^m^ ^ Uj, 



(4.1) 



where lij is defined in 2.3, 2.4. The benchmark point gives mj, = 0.05 eV, consistent with 
experimental limits on the sum of neutrino masses ~ 1 eV. Furthermore, varying the Dirac 
Wino mass between 500 GeV and 1 TeV does not change the statistics from the neutrino 
mass and mixing cuts. 

Experimentally, the observable sc which corresponds to the ^13 angle in the neutrino 
mixing matrix U is the least constrained [34] . The randomly generated sample of neutrino 
mass matrices tends to have a sc distribution that peaks towards zero mixing, as shown 
in Fig. 2. A non-zero measurement of ^13 at Daya Bay or at another similar experiment, 
for instance, would support an anarchic scenario, but even if ^13 were constrained to be 
smaller than 0.01, that would not exclude it. 




0.01 

a 

0.00 




sin^e. 



sterile 



■ 2 „ 

sm 9. 



sterile 



Figure 3: Neutrino mass cuts make the sneutrino LSPs more sterile: (i) Run I (ii) Run II. The 
black plot satisfies the cuts and the light green curve is a random sample that docs not satisfy the 
neutrino constraints. 



One very interesting feature of our scenario is that the neutrino mass constraints shed 
light on the active/sterile mixing of the LSP sneutrino. In particular, the sneutrino LSP 
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Figure 4: e-fi-r mixing in the active sector for the sneutrino LSPs that satisfy neutrino mass cuts: 
(i) Run I (ii) Run II 



tends to be more sterile, as seen in Fig. 3. This bias toward a mostly sterile LSP is mostly 
due to the R cut. Thus, we do find that at least the requirement of a mass hierarchy in 
the neutrino sector seems to prefer a sterile-dominant LSP, which is in accordance with 
expectations from relic abundance and Z-width constraints. Furthermore, in the active 
flavor sector, the sneutrino LSP tends to be predominantly ^ or r flavored as seen in Fig. 
4. 

4.2 Relic abundance 

The scatter plots of Fig. 5 show the Oda/^^ values for the random matrices of Runs I 
and II. The funnel shapes at the Z and Higgs poles are as expected. Beyond this, the 
most striking features of these plots is that the relic abundance constraint selects a very 
"ordinary" region of parameter space. We needn't go into the Higgs or Z poles or out to 
the edges of the parameter scan to find LSPs with the appropriate relic abundance. 

As discussed earlier. Table 1 shows that a random set of sneutrino mass matrices 
is as likely to satisfy the relic abundance constraint as one satisfying the neutrino mass 
constraints. The Q^oKih'^ constraint is furthermore essentially flavor-blind and so does little 
to change our expectations for the flavor content of the LSP sneutrino, compared with what 
we found after the neutrino cuts alone. This can be seen by comparing Figures 4 and 6. 

The distributions of LSP masses and mixing angles for the sets of mass matrices that 
satisfy both the neutrino mass and relic abundance constraints are plotted in Figures 7 and 
8. The most important result is that the LSP typically has very little electron component, 
and is dominantly a mixture of and r. This means that it is likely that cascades to the 
LSP will be /i and r rich, motivating the collider analyses in section 5. 
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sin ee/siu Sactk 



sin ^e/siu ^active 




Figure 6: e-^-T mixing in the active sector for the sneutrino LSPs that satisfy neutrino mass and 
rehc abundance constraints: (i) Run I (ii) Run II 



4.3 Direct detection 

XENONIO [41] and CDMS [40] put the strongest constraints on spin-independent WIMP- 
nucleon scattering. Because the Z-exchange contribution to the scattering amphtude can 
be strongly suppressed in the presence of lepton-number violation, we require that the 
Higgs-mediated contribution satisfy the bounds from these experiments. Fig. 9 shows the 
distributions of LSP-nucleon cross sections for the sets of mass matrices that satisfy the 
neutrino, relic abundance, and invisible Z-width constraints. 

Using the same sets of matrices, distributions of sneutrino mass ratios are shown in 
Figures 10 and 11. The vertical gaps in the plots are due to the fact that sneutrino LSPs 
near the Z or Higgs poles rapidly annihilate to leave a very low relic abundance. These 
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Figure 7: Mass and mixing angle distributions for the 343 sneutrino LSPs from the sample that 
satisfies neutrino mass, relic abundance and invisible Z width cuts in Run I. The LSP mass is in 
GeV. 
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Figure 8: Mass and mixing angle distributions for the 144 sneutrino LSPs from the sample that 
satisfies neutrino mass, relic abundance and invisible Z width cuts in Run II. The LSP mass is in 
GeV. 



plots distinguish between the sneutrino LSPs that do and do not satisfy direct detection 
bounds. While a significant fraction do in fact satisfy the current direct detection bounds, 
the next generation of direct detection experiments will probe down to the 10~^^cm^ level, 
and can therefore exclude much of the remaining parameter space. 

5. Collider phenomenology 

The LHC phenomenology of mixed sneutrinos in the MRSSM is potentially very rich. If 
we expect large flavor mixing in the sneutrino mass matrix, we should also expect left- 
handed slepton mass matrices with large off-diagonal entries. We will see that cascades 
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'"lsp (GeV) m^sp (GeV) 

Figure 9: Distributions of I'l-nucleon spin-independent cross sections for Higgs-mediated scatter- 
ing, for the samples that satisfy neutrino mass, reUc abundance and invisible Z width cuts. The 
regions above the thick (blue) curve is excluded by CDMS and XENONIO. In Runs I (left) and II 
(right), 69 and 12 LSPs satisfy the direct detection constraints, respectively. 



involving left-handed sleptons can consequently produce opposite-sign dileptons that reveal 
this lepton flavor violation. 

Because experimental constraints on flavor violation are the strongest for the first two 
generations, e-[i flavor violation associated with sleptons produced at the LHC would be 
a particularly striking signature of this scenario. Moreover, as we have learned in the 
previous section, this model tends to give produce an LSP with a significant /i-component, 
motivating inclusive signatures of this feature. In this section we show that in certain 
regions of MRSSM parameter space, these signatures will indeed be observable. Recent 
work has shown detectability of a r flavor violating signal [46], and other related work on 
flavor violation includes [14, 47, 48]. 

5.1 Slepton mass matrix 

The slepton mass matrix is obtained from the highly mixed sneutrino mass matrix, by en- 
forcing the constant D-term mass difference in the left-handed sector. Upon diagonlization, 
the left-handed mass eigenstates are related to the gauge eigenstates by li = Uiele where 
Uie diagonalizes the mass matrix, i = 1,2,3. Unlike in the MSSM, where the left-handed 
selectron and smuon are nearly degenerate, the left-handed sleptons have a mass hierarchy 
that can lead to an observable flavor violating signal. 

Although the right-handed slepton mass matrix may also be highly mixed, we take it 
to be diagonal for simplicity's sake. The e-fi flavor-violating signal we study below can 
arise just as easily in the more general case with mixed left- and right-handed slepton mass 
matrices. 

5.2 Neutralino and chargino masses 

The MRSSM has 4 neutralinos in addition to the MSSM neutralinos. In the absence of a 
jj, term, the and pair up with iJ^ and iJ^ respectively, HuH^Ru and ^dH^R^^. Fur- 
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Figure 10: Sneutrino mass ratios for the 343 sets of matrices that satisfy neutrino mass, rehc 
abundance and invisible Z width cuts in Run I. The grey circles and black asterisks represent 
sneutrino LSPs that do not/do satisfy XENONIO and CDMS bounds. 



thermore, the bino and wino get Dirac masses with two additional neutrahno fields. These 
Dirac gauginos induce radiative contributions to the sfermion masses-squared that are a 
factor ~ g(^/(167r^) smaller than the gaugino masses-squared themselves. The charginos 
are similarly paired with the charged components of the additional fields. 

5.3 e-fi flavor violation 

A decay chain involving left-handed sleptons and ending in vi can provide a signal of lepton 
flavor violation. An important example for the sample parameter point we study below is 
the chain x? hi vivl, where xi is the lightest neutralino and li is the lightest left- 
handed slepton (see figure 12). As the VF-boson in this decay chain has flavor-universal 
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Figure 11: Sneutrino mass ratios for the 144 sets of matrices that satisfy neutrino mass, relic 
abundance and invisible Z width cuts in Run II. The grey circles and black asterisks represent 
sneutrino LSPs that do not/do satisfy XENONIO and CDMS bounds. 



couplings, flavor violation occurs at only the first vertex, Xi — — I- The sneutrino vertex 
itself does not contribute to the flavor- violating eflFect. The signal we will focus on therefore 
provides information on the flavor structure of the lightest left-handed slepton, li. 

Having a larger wino than bino component in the lightest neutralino/chargino increases 
the coupling to the left-handed sleptons and favors more flavor violating events, given that 
we are taking the right-handed slepton mass matrix to be flavor-diagonal. The e-/x flavor- 
violating signal is also enhanced if the lightest neutralino couples to the leptons dominantly 
via gauge couplings rather than Yukawa couplings, so that decays to r's are suppressed. 
So, it is advantageous if x5 has a larger than component. These features are realized 
in the parameter point that we adopt for our detailed analysis. 
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Figure 12: A decay chain with potential lepton flavor violation 



5.4 Particle spectrum 

We now describe the most important details of our sample parameter point. The parame- 
ters that determine the neutralino and chargino properties are: 

tan /3 = 10 fj,u = 348 GeV /x^ = 359 GeV 
mDwino = 1 TeV mDuno = 2 TeV 

where mD^ino and mDhino are the Dirac masses of the wino and bino. 

The bino is heavier than the wino in this setup, but this is not at odds with gauge 
unification. The new fields in Dirac gaugino models pair up with additional "bachelor" 
fields to form complete multiplets and perturbative unification is possible for gauge groups 
SU{3)^ [4] or SU{5). The relative sizes of the gaugino masses are not entirely determined 
and depend on the interactions with the "bachelor" fields. Thus a bino heavier than the 
wino is still consistent with unification of gauge couplings. 

We adopt one of the randomly generated sets of sneutrino mass matrices that passes 
the relic abundance, neutrino mass, Z-width and direct detection cuts described in Section 
4. The left-handed slepton mass matrix is then fixed accordingly. The sparticle masses 
and the flavor structure of the LSP sneutrino are given in Table 5. In the parameter 
point considered, the lightest neutralinos and charginos are predominantly of Higgsino- 
like and their masses are much lower than the TeV scale. However, these light masses are 
consistent with the anarchic sfermion mass matrices since flavor observables in the MRSSM 
are suppressed by the gaugino masses, which are ©(TeV). 

We include only the lightest four neutralinos, the lightest two charginos, and the 
lightest three sneutrinos in our analysis as the heavier ones are not produced appreciably. 
Truncated in this way, the spectrum is easy to compare with an MSSM spectrum. Note 
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Sparticle 


Mass (GeV) 




1 TeV 


m-g 


1 TeV 


Xl,2 


347 


A3, 4 


359 


m-+ 


346 


xt 


359 




ZOO 








311 




375 




78 




124 




227 



v\ flavor structure 


Sin^^active 


0.11 


Sin^^stcrilc 


0.89 




-0.05 




-0.324 




0.084 



Table 5: Superpartner masses and sneutrino LSP flavor structure for the chosen parameter point. 
V relates the gauge and mass eigenstates. 



that the neutralinos and charginos have almost degenerate masses and e.g. decays from a 
heavier to a lighter neutralino are not allowed kinematically. This feature is not essential 
for the e-[x flavor-violating signal, however. 

The masses of the lightest neutralinos and charginos are approximately set by the 
values of \Xu and /i^. These parameters are constrained by electroweak symmetry breaking 
and cannot be too large. Increasing /^^ relative to [id increases the number of r decays due 
to the Yukawa coupling of if^, the lightest neutralino. Increasing the masses of the lightest 
neutralinos and charginos opens up new decay channels through the other mixed sleptons 
(rather that just through /i), resulting in a richer LFV signal with multiple sources of flavor 
violation. 

For this parameter point the LSP v\ is primarily sterile with a large /i component in 
the active flavor sector. More important for the collider phenomenology is the fact that the 
lightest left-handed slepton also has a large [i component and tiny e component. Decays 
to left-handed sleptons primarily go through l\ since it is the lightest, so we should expect 
an excess of muons over electrons. On the other hand, decays to right-handed sleptons 
proceed via f2(= = Ir^) due to the large Yukawa coupling. 

The mixing angles for the MSSM components of the lightest four neutralinos are given 
in Table 5.4. The lightest neutralino Xi 2 mainly and has a larger wino than bino 
component. The wino component can be increased by increasing the Dirac bino mass. 
As already mentioned, a larger wino component increases the coupling of neutralinos to 
left-handed sleptons to produce a stronger LFV signal. 
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Neutralino 






m 


u 


Al,2 


T 0.00275 


± 0.02212 


0.0079 


-0.7041 


vO 
A3,4 


± 0.00025 


=F 0.00204 


-0.707 


-0.0085 



Table 6: The composition of the four Hghtest neutrahnos for the chosen parameter point. The 
mass eigenstates also contain mixtures of the Dirac partners of 5°, , H'^, and 



5.5 Branching Ratios 

The relevant branching ratios for the chosen parameter point are given in Table 7. We 
will be particularly interested in flavor violation in events with ee, /i/i, and e/i pairs, which 
arises dominantly from the 2 decay chain depicted in Figure 12. 



X"2 ^ ll,T2 

^0 0.246/0.14 7 , , . 0.66 

Xi,2 > hif^/r) > Wui 



0.11 



> Dili/ 



_0 0.005 - 0.495 - 
Xl,2 > > V2TV 



0.425 r hdccays 
> l\TT 



X3,4 ^ ll,T2 



.0 10-^0.423 7 , I . 0.66 
X3,4 > h^^lr) > Wv 



0.11 , 
> V\lV 



~0 0.56 ~ 0.495 ~ 
X3 4 > T2T > V2TV 



0.425 r iidocays 
> l\TT > 



Xi decays 



-+ 0.935 r 
Xl > 



X2 decays 


0.56 f2decays 

X2 > > ... 

3.6x10-4/0.215 . , , - 0.1 , 
X2 ^ J^2(/^/t) ^ + « 


1.8xl0-*/0.014 _ , , , 

X2 

1.8xl0-*/0.215 - , , s 0.036 , , 

X2 ^ J^slAt/Tj ^ z^i/+/ 



Table 7: Branching ratios for the neutralino and chargino show flavour violation at the left-handed 
slepton, Zi, vertex 



As shown in Table 7, Xi 2 decays to and Zir with about the same branching ratios. 
In the absence of Yukawa couplings, the branching ratio to would be much larger 
because l\ is dominantly /u-flavored. However, the r Yukawa coupling makes tau production 
competitive with muon production, because x? 2 has a small, non-negligible ff^ component. 
In contrast, x-j 4 has a dominant component, so X3 4 decays are driven by the r Yukawa 
coupling, leading to far more taus than muons. 

As for the right-handed sleptons, the decay f2 — )• Iitt via off-shell neutralinos is 
dominant over decays to lier or lifir. The large r Yukawa coupling means that f2 couples 
more strongly to X34 ~ than to X12 ~ ^u- ^o, the three-body decay is dominantly 
mediated by X34) which preferentially produces r's. For the parameter point chosen, 
though, f2 decays are not important for the analysis we describe below. 
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5.6 Analysis 

To study the flavor- violating signals of this parameter point at the LHC, we calculate 
branching ratios for 2- and 3- body decays in CalcHEP (with modified vertices), and then 
pass these tables to Pythia 6.4 [49], which we use to simulate full SUSY production and 
parton showering. We also use Pythia to generate a sample of ti background events. 

We generate 96,020 SUSY events, corresponding to 100 fb^^ of integrated luminosity. 
Because the gluinos are Dirac fermions, we turn off LL, RR, LR* squark-squark produc- 
tion. The SUSY production modes g-qL,R, Ql-Qr, g-g, and q-q* account for 63.3%, 9.9%, 
9.3% and 7.1% of the total SUSY events generated. The remaining 10.3% comes from 
neutralino/chargino production. 

To suppress the background from standard model processes, we keep only those events 
with 

• at least two leptons with pT > W GeV and \r]\ < 2.4 

• Ept > 1500 GeV, 

where X^pt includes jets with pT > 20 GeV, leptons with px > 10 GeV and |?/| < 2.4, 
and photons with pT > 10 GeV. After these cuts, the leading SM background is tt. The 
number of tt events that pass these cuts after 100 fb~^ of integrated luminosity is 997, 
versus 20, 173 events from SUSY production. 

The analysis for e-fi flavour violation is not qualitatively affected if the dominant 
production mode, g-qL,R, becomes unavailable, say in the case of a very heavy gluino. For 
a 2 TeV gluino, the number of SUSY events generated is 20, 650 and the number of SUSY 
events that pass all cuts is 2, 817. Although the signal is scaled down, the flavor-subtraction 
techniques employed efficiently reduce the tt background and a flavor- violating signal is still 
observable. 

The distributions of dilepton invariant masses have various associated endpoints. The 
kinematic endpoint for the dileptons in the type of decay chain shown in Figure 12 {x^ — )• 
III — )• Dul) is 

mf — m?- + mlj + \/ {mf — m?- + )2 — (2mjmwY 
— ^ =. (5.1) 

TTl? — m? + 771,2, _ . (m- — 777,2 _|_ 77j2^)2 _ (2777,^-777,^)2 

The kinematic endpoint for the taus in the decay chain — )• — )• vvt has a simpler 
expression, as it involves one two-body decay and one-three body decay, rather than three 
two-body decays: 

/ (777?,! — 777? 1(777? — 777?) 
^raa. ^ ^^^^ ^ (5_2) 

Numerical values of endpoints associated with cascades involving li and f2 are given 
in Table 8. 



777, 
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Decay Channel 




Xi,2 ^ h ^1 


257 


X3,4 h Z>1 


269 


Xl,2 T2 Z>2 


207 


X3,4 T2 Z>2 


222 



Table 8: Kinematic endpoints for decays via h and f2. 

We will focus on ee fi^ and e/i opposite-sign dilepton pairs, which singles out the 
Xi,2 h ^ i^i decay chain (x3,4 — )• — )• t'l occurs much less often). Taus produced in 
f2 chains could potentially contaminate the e-/i signal that interests us when they decay 
leptonically, but these decays are flavor universal and are effectively subtracted out when 
we consider distributions that are sensitive to e-fi flavor violation. In our analysis, we 
include the effect of the decays of rs that are produced in the SUSY decay chains and ti 
background. 

A number of strategies can be applied to gather evidence for e-/i violation. The prob- 
abilities of having an electron or a muon come from the first vertex in the decay chain of 
Figure 12 are directly proportional, respectively, to the e and /x contents of li. Because 
electrons and muons are produced with equal probability by the W^-boson in the chain, 
only if the e-/i mixing in /i is 1 : 1 will the number of ee and fifj, events be equal on average. 
By considering the ratios of the numbers of ee, e/i, and fifi events, information about the 
flavor violating slepton vertex can be extracted. 

Event counts for opposite-sign ee, e/i, and /i/z pairs with invariant masses in the range 
180 GeV < mil ^ 290 GeV are given in Table 9. These counts suggest that the lightest 
left-handed slepton has a larger fi than e component. The mass range is chosen to ensure 
that the events arise predominantly from h decays. 



Distribution 


Multiplicity 


ee 


216 


mi 


1,711 


eji 


1,601 



Table 9: The number of opposite-sign ee, /t/t and e/t events in the mass range 180 GeV < mil < 
290 GeV. 

Subtracting ee{OS — SS) invariant mass distribution from the ^^{OS — SS) invariant 
mass distribution (Fig. 13) not only shows that there are more /i than e events, but also 
reveals a clear endpoint near the expected value. The tt background is suppressed by the 
flavor subtraction and the SS subtraction reduces combinatorial background in the SUSY 
events. The derived endpoint is calculated by a 2-line fit on the data, as depicted by the 
straight lines in the plot. 

In the absence of fiavor violation, there should be twice as many efi{OS) events as 
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Figure 13: An excess of /i/x over ee events indicates e-fj, flavor violation 
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Figure 14: The ratio of jifi to efj, events is not 1:2, suggesting e-/i flavor violation 



fj,^fj-^ events. An excess after subtracting e// events from twice the number of /i/i events 
is therefore another indicator e-^ flavor violation. In Table 9 we see that the numbers of 
fJ^~^fJ'~ and efi{OS) events are almost the same, rather than different by a factor of two, 
and the plot in Fig. 14 shows the relevant invariant mass distribution. The SS subtraction 
again reduces combinatorial background in the /U and e^{OS) SUSY events. 

The e-/i flavor violation can also be seen in the e/i events alone. In the decay chain 
of Figure 12, the lepton from the neutralino-slepton vertex typically has a higher pT than 
the lepton arising from W decay. The flavor distribution of the harder lepton can thus be 
used to learn about the flavor content of h. 
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Figure 15: distributions for events with highest px e subtracted from highest px fJ^ events. 



From the total sample of 20,173 events that pass all px cuts and contain at least two 
leptons, the numbers of events whose highest-pr lepton is of e or type are given in Table 
10. Since this is an inclusive study of all possible SUSY decay chains, exact mixing angles 
cannot be derived from this analysis. However, the numbers clearly indicate that hard 
muons are favored over hard electrons. 



Highest pt lepton 


Multiplicity 


e 


5,391 




14,782 



Table 10: The numbers of events passing cuts with the hardest lepton being e or /z. 



By subtracting the dilepton invariant mass distribution for efi events in which e is the 
hardest lepton from the distribution for e/i events in which fi is the hardest lepton (Fig. 
15), a clear event excess indicating e-fi flavor violation can be observed. 

That the derived endpoint does not match the endpoints seen in the previous plots 
is not surprising given that the bias toward harder muons does not persist throughout all 
regions of phase space. For example, if the sequence of two-body decays depicted in Figure 
12 proceeds in a way that maximizes the invariant mass of the two charged leptons, then 
for the parameters we've adopted it turns out that the lepton from the first vertex (the one 
more likely to be a muon than an electron) is in fact slightly less energetic that the lepton 
from the second vertex, in the rest frame of the neutralino. 

A more sophisticated analysis might improve discrimination of signal and background 
and the readability of kinematic endpoints. It would also be interesting to employ some 
level of detector simulation to explore how the types of analyses we have described might 
work in a real detector environment. In particular, the different efficiencies for detecting 
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electrons and muons may complicate the flavor subtraction technique we have used to 
isolate the flavor-violating signal. 

6. Conclusions 

Dark matter remains one of the most exciting topics in particle physics, and one in which 
some of the most significant experimental progress will be made over the next few years. 
We have revisited the topic of mixed sneutrino dark matter in the context of the Minimal 
R-symmetric Supersymmetric Standard Model (MRSSM), focusing in particular on the 
connections between neutrino masses and mixings, dark matter, and collider phenomenol- 
ogy in this framework. 

Because flavor violation is expected to be large in the MRSSM, the anarchic generation 
of large neutrino mixing angles flnds a natural home in this theory. The small Majorana 
gaugino masses allow the radiatively generated neutrino masses to be of the appropriate 
size, while still allowing enough lepton-number violation for the LSP sneutrino to evade 
strong direct detection constraints on scattering via Z-exchange. We find that the region 
of neutrino parameter space defined by the current experimental constraints is consistent 
with the idea that the neutrino mass matrix is generated by an essentially random set of 
sneutrino mass matrices. 

Data from neutrino experiments give some insight into the expected properties of the 
LSP sneutrino. First, the LSP tends to be dominantly sterile. This is to the benefit of the 
model, as this same region of parameter space is preferred by the measured relic abundance. 
Second, the LSP tends to have significant /x and r components. The consequence of a large 
fx component in the LSP is that studies of fxfj,, ee, and e/x invariant-mass distributions 
can give conclusive evidence of lepton flavor violation at the LHC. In particular, inclusive 
signatures of opposite sign dileptons can reveal strong flavor violation that would not 
appear in the MSSM, even with small non-degeneracies in the first two generations. 

Direct detection experiments have a great possibility of probing this model. While 
the model provides a natural basis for inelastic dark matter, it also predicts a sizeable 
elastic scattering cross section 10^^'^cm^) and should give detectable signals in the 
next generation of dark matter direct detection experiments. Should such a scenario be 
discovered, we may be able to probe the origin of the mysterious properties of the neutrino 
mass matrix at the weak scale. 
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